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DESIGN EQUATIONS OF SWEPT FREQUENCY SPECTRUM ANALYZERS 
FOR IN-FLIGHT VIBRATION ANALYSIS 
by Edward F. Miller and David L. Wright 
Lewis Research Center 

SUMMARY 

The equations for the design of swept frequency spectrum analyzers are presented 
with the intended application being in-flight analysis of vibration data. Primary emphasis 
is given to developing the optimum combination of the frequency- time function and the 
analyzing filter bandwidth. The main criterion used is minimum sweep time coupled with 
acceptable resolution and accuracy. Constant, linear, and second-order sweep rates are 
discussed, and the frequency-time equation is developed for each of these three types. 

The analyzer using a second- order sweep rate is shown to have the shortest sweep time, 
but it requires an expanding bandwidth filter and a decreasing averaging- time power de- 
tector. The linear sweep rate analyzer is shown to require less bandwidth per channel 
and is much less complex. The constant sweep rate analyzer is shown to have either very 
long sweep times or poor frequency resolution; however, all three types of analyzers 
have bandwidth savings over real-time transmission of wide-band data. 


INTRODUCTION 

A continuing problem in communication system design is transmitting many channels 
of wide-band data over a limited bandwidth transmission system. The problem is most 
critical in transmission from a spacecraft or a launch vehicle, where increased band- 
width requires more weight and power, and consequently, reduced payload. The partic- 
ular problem that led to this investigation was the need to transmit from launch vehicles 
the results of a large number of vibration measurements. 

Currently, NASA uses primarily the IRIG FM/FM (Inter-Range Instrumentation 
Group) standard telemetry format for the transmission of analog vibration signals from 
vehicles to the ground, and an extensive network of ground stations has been built to re- 
ceive data in this format. Often, the number of vibration measurements desired exceeds 



the capacity of the IRIG FM/FM telemetry system. There are several ways to solve this 
problem. One method is to change the transmission system. Two transmission systems 
capable of transmitting many channels of real-time, wide-band data from space vehicles 
are DSB/FM (double-sideband/frequency modulation (ref. 1)) and SSB/FM (single- 
sideband/frequency modulation (ref. 2)). Another alternative is to use some form of data 
reduction before the data is transmitted. On-board data reduction reduces the transmis- 
sion bandwidth required and still allows the use of the IRIG FM/FM transmission system 
for the resulting low bandwidth channels. Spectrum analysis is the method of pretrans- 
mission data reduction that is considered in this report. 

Spectrum analysis of a signal often gives sufficient information about the real-time 
signal. In many cases, the final form of vibration data is a plot of the power density of 
the signal as a function of frequency (or power spectral density). This form of the re- 
sult is important because it can show the frequencies present in the vibration source as 
well as the amplitudes and frequencies of the dominant vibration modes of the vehicle 
structure. However, the process of spectrum analysis is irreversible since no phase 
relations between the frequency components of a given signal are obtained; thus, spectrum 
analysis loses the actual time history of the vibration. In addition, relative phase meas- 
urements between two time functions cannot be made by comparing their spectrums. Even 
with these limitations spectrum analysis remains useful. 

Theoretical work in spectrum analysis (refs. 3 and 4) was followed by the consider- 
ation of swept frequency spectrum analyzers with a constant sweep rate and a constant 
bandwidth filter (refs. 5 and 6). Commercially available airborne spectrum analyzers 
are of this type; however, they have the drawbacks of poor resolution and/or long sweep 
times. For random data, a spectrum analyzer with a linear sweep rate has better reso- 
lution and/or shorter sweep time than the constant sweep rate analyzer. The linear sweep 
rate analyzers have been used (to the authors* knowledge) only as ground equipment. A 
spectrum analyzer with a sweep rate proportional to frequency squared was mentioned by 
Bendat (ref. 7), but the frequency-time function was not derived. 

In this report, the basic theoretical background for spectrum analysis is presented. 
The frequency-time functions are derived for the various types of swept frequency spec- 
trum analyzers. The development progresses from the simple case of a constant sweep 
rate analyzer operating on periodic signals through to the progressively more complex 
operations of constant, linear, and second-order sweep rate systems analyzing random 
signals. For each type of frequency sweep, the attainable performance, as given by the 
frequency-time function, is compared with the following desired specifications for air- 
borne spectrum analyzers (the list was compiled from discussions with those who use and 
reduce vibration data): 
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Frequency range near zero to several kilohertz 

Allowable frequency error ±10 percent of the frequency being analyzed 

Minimum resolvable bandwidth ±10 percent of the frequency being analyzed 

Allowable spectrum amplitude error ±10 percent 

Sweep time of the order of seconds 

These desired specifications are shown to be incompatible and cannot be obtained in 
the analysis of random data. Even periodic data can only be analyzed with less than the 
desired accuracy in such a short sweep time. 

In those cases where long sweep times or reduced accuracies are tolerable, the in- 
flight spectrum analyzer may be used. The bandwidth savings resulting from the use of 
spectrum analyzers in these cases is discussed in the concluding section of this report. 

A limitation of the swept frequency spectrum analyzers discussed in this report is 
that it is not practical to extend their lower frequency limit to ’’near zero, '* as listed in 
the desired specifications previously given. The frequency-time equations derived in this 
report show that low-frequency spectrum analyzers have extremely long sweep times and 
require very narrow filters to obtain good frequency resolution. The swept frequency 
spectrum analyzers have a practical low-frequency limit on the order of 100 hertz. The 
low-frequency part of the signal is better handled by separate transmission of that portion 
of the time-varying signal followed by spectrum analysis on the ground. 


BACKGROUND TO PROBLEM 

Before proceeding to the details of the analysis, some necessary background infor- 
mation will be discussed. 


General Background 

Characteristics of expected signals. - Since this report is concerned with the specific 
problem of transmitting vibration data, the assumed characteristics of such data will be 
discussed. First, the source causing the vibration is assumed to be a random noise 
source. The mechanical system coupled to this source acts as a filter that shapes the 
frequency content of the resulting vibration. The mechanical Q, or sharpness of reso- 
nance, is assumed to be less than 10; thus, if the spectral peaks are to be resolved, 
analyzing filters with values of electrical Q much greater than 10 would be required. 
Moreover, the maximum mechanical Q is assumed to be independent of the resonant 
frequency. Thus, the bandwidths of the resonances become progressively wider as the 
frequency increases. 


3 



In order to obtain meaningful data, the signal being analyzed must be stationary 
during the time of observation; that is, the probability distribution of the signal level 
must remain constant. The vibration signal will remain stationary as long as the source 
causing the vibration remains statistically invariant. Examples of events causing nonsta- 
tionarity would be turning on or turning off a rocket engine, changes in thrust, or 
changes in air turbulence surrounding the vehicle. 

For the swept frequency spectrum analyzer, the minimum period of stationarity of 
the signal limits the total sweep time. When the signal is nonstationary, the spectrum 
changes with time. Thus, for the swept frequency spectrum analyzer the signal must be 
stationary for times on the order of the total sweep time. The minimum period of sta- 
tionarity can be somewhat less than the total sweep time since interpolation between suc- 
ceeding spectrums can give estimates of the spectrum between sweeps. 

Types of spectrum analyzers . - There are two common types of analog spectrum 
analyzers. The simplest method of spectrum analysis is to measure the power from each 
filter in a group of band-pass filters whose center frequencies are spaced so that the filter 
pass bands cover the frequency range of the input signal. The power output of each filter 
is an approximation to the average power spectral density of the signal within the filter 
pass band. The swept frequency spectrum analyzer is the second common type of spec- 
trum analyzer and is the subject of this report. Its operation will be discussed in de- 
tail. 


The analyzer shown in figure 1 is one of several ways of implementing the swept fre- 

Basically, the analyzer in 


quency analyzer. The frequency diagram is shown in figure 2. 
figure 1 translates the signal spectrum to a higher frequency using suppressed carrier 
amplitude modulation. At the high frequency, a single band-pass filter selects a narrow 
band of frequencies corresponding to a similar narrow band of frequencies of the input 

spectrum. A sample of the power out- 

Spec- 
trum 
output 



Figure 1. - Block diagram of typical swept frequency spectrum analyzer. 



Figure 2. - Frequency diagram of typical swept frequency spectrum analyzer. 
Lowest frequency to be analyzed f L ; highest frequency to be analyzed f H ; 
mixing frequency f s ; center frequency f c . 


put of the filter is an approximation to 
the power spectral density within the 
corresponding signal frequency incre- 
ment. The filter is effectively swept 
through the signal spectrum by sweep- 
ing the mixing frequency f . A disad- 
vantage of this analyzer is that both the 
filter and the power detector must have 
time to change their output values as 
they are swept through the spectrum. 

A major objective of this report is to 
show that the proper choice of the time 
dependence of the frequency sweep will 
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optimize the sweep time and the resolution for. the swept frequency spectrum analyzer. 
Before deriving the frequency-time equations for the constant, the linear, and the second- 
order sweep rate spectrum analyzers, some of the basic theory of spectrum analysis will 
be presented. 


Theoretical Background 

The basic restrictions on swept frequency spectrum analysis will be shown by first 
considering the spectrum analysis of periodic signals. Then the more complex case of 
analysis of random signals will be considered. 

Spectrum analys is of periodic signals . - Periodic signals are representable by a 
Fourier series of harmonically related sinusoidal waves. The basis for such representa- 
tion is well known and will not be discussed here (see ref. 8 for a discussion of Fourier 
series). 

The spectrum for a periodic signal has discrete components whose amplitudes are 
constant with time. Such a spectrum may be resolved by a swept frequency spectrum 
analyzer with a sharp band-pass filter whose bandwidth is less than the minimum spacing 
between frequency components of the spectrum. Thus, 

B f <-1 (1) 

1 T„ 


where 


Bj effective filter bandwidth (bandwidth of the ideal rectangular filter that passes the 
same power as the real filter when both are excited by white noise) 

Tp period of the signal 

Furthermore, the sweep rate of the analyzer must be slow enough to allow the filter 
response to come up to almost full value before the analyzer moves on to another fre- 
quency. A band-pass filter excited by a sinusoidal signal at its center frequency has an 
output response envelope with a time constant given by 



( 2 ) 


where T„ is the time constant of response. 
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Also, the output of the filter must be measured by a detector whose averaging time is 
long compared with the period of the filter output, which is a sine wave at the filter center 
frequency. Thus, the averaging time is given by 

T d »J- (3) 


where 

detector averaging time 
f filter center frequency 

Using equations (1), (2), and (3), one can design a spectrum analyzer for periodic signals. 

Spectrum analysis of random signals . - The measurement of the spectrum of a ran- 
dom signal is not as straightforward. By its nature, a random signal has no periodicity, 
its value as a function of time cannot be precisely predicted, and, at best, only a prob- 
ability distribution for its values can be specified. Analyzing such functions requires 
some of the techniques of generalized harmonic analysis and probability theory. A brief 
introduction to the required techniques follows (for a thorough discussion see refs. 9 to 
11 .). 

The power spectral density (or mean square spectral density) of a stationary random 
function of time is defined as the Fourier transform of the autocorrelation function R(r). 
Let x(t) be the random function. Then 


R(t) = lim 

T— °° 2T 



x(t) x(t+r)dt 


and the power spectral density S(f) is given by 


(4) 


S(f) = f R(r)e _;i27rfT dT 

*/_ on 


(5) 


(Symbols are defined in the appendix. ) The function S(f) has the property that 


L 


S(f)df = Mean square value of x(t) = 

00 T 


lim l_ r 

T ~°° T J Q 


x 2 (t)dt 


( 6 ) 


Also, 
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( 7 ) 



Sffldf . P ab 


where P ab is the average power in x(t) contained between frequencies f and f^, and 
between frequencies -f and -f^, where f^ > f & > 0. If S(f) is nearly constant between 
f and f b , for example S ab (f), then S ab (f) may be approximated by 



(f) 


ab 


2 < £ b - y 


( 8 ) 


where the caret above the S indicates that this is an estimated value. 

The spectral density previously defined is two sided; that is, S(f) is nonzero for both 
positive and negative frequencies. In engineering applications, spectral densities are 
usually considered one sided (nonzero only for positive frequencies). Hereinafter, the 
spectral densities will be considered one sided. Equation (8) becomes 


ib< f > 



m 


This equation is the basis of the most common method of spectral analysis, the use 
of a narrow band-pass filter to measure the power in x(t) between the frequencies f 
and f^. If the effective filter bandwidth is B^, then 


s ab< f > 


ab 


( 10 ) 


where P ab is the average power out of the filter whose effective bandwidth is B^; that 
is, 



lim 

T— 



Xj (t)dt 


(U) 


where x^(t) is the output of the filter. 

Rewriting equation (10), incorporating equation (11), and considering the limiting 
case yield the following expression for the exact spectral density: 
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I 


( 12 ) 


S «c> ' 


lim 

B { -*0 


lim 



where 

f center frequency of the filter 

V 

S(f c ) exact value of the power spectral density at frequency f c 

For practical reasons, Bj cannot go to zero; however, a small value of B^ will 
give a close approximation as long as T is sufficiently large. 

The estimated value of S(f ) can be written as 

s(f c )=-i_ r Xf 2 (t)dt d3) 

B f T Jq 

In words, S(f„) is the estimate of the average power spectral density in a bandwidth B f 
centered about f„. S(f r ) is an unbiased estimator of the spectral density at f only if 
the spectrum is constant in the bandwidth Bj. (However, for nonflat spectrums, S(f c ) is 
an unbiased estimator of the average spectral density within the bandwidth B^. ) Since 
S(f ) is only an estimate of the true value of S(f ), the measured values of S(f ) will have 
statistical variation about their average value. 

When equation (13) is implemented, T is the averaging time of the power detector on 
the output of the filter. If the power detector implements equation (11) for finite T 
through an integrate and reset mechanism, true averaging is performed and T is the 
time between integrator resets. However, if a resistance-capacitance averaging circuit 
is used, T is two times the RC time constant (ref. 7, section 4). 

The probability distribution of S(f ) will indicate the expected range of values of these 

v A 

estimates of the spectral density. To determine the probability distribution of S(f c ) con- 
sider the following: 

(1) From sampling theory, it can be established that the maximum number of inde- 
pendent samples obtainable from a band-limited white noise signal of bandwidth B f and 
duration T is n = 2B^T (ref. 12). 

(2) The variance estimate of a normal distribution with zero mean value is a function 
of the chi-square random variable with N-l degrees of freedom, where N is the num- 
ber of samples taken to determine the sample variance (ref. 13). Specifically, 

^(N-l)=x|.i < 14 > 

a 
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where 

sample variance 
true variance 

chi-square random variable with N-l degrees of freedom (the chi-square distri- 
bution is tabulated in statistical tables) 

(3) For certain types of random functions, those which satisfy the quasi-ergodic hy- 
pothesis (ref. 11, ch. 7), the time average of the function equals the ensemble average; 
that is, 


2 

X N -1 


lim 


T— °o 



x 2 (t)dt = 



p[x(t)]x 2 (t)dx 


(15) 


where p[x(t)] is the probability distribution of x(t). Equation (15) states that the mean 
square value of x(t) equals the second moment of x(t). The random functions considered 
will be assumed to satisfy the quasi-ergodic hypothesis. 

(4) For a random variable x(t), 

E[x 2 (t)] = Var x(t) + |E[x(t)]| 2 (16) 

f A \ r A i ( n 

Second Variance Mean value 

moment of x(t) of x(t), quan- 
of x(t) tity squared 


Applying equation (6) and steps (1), (2), (3), and (4) leads to the following conclusion: 
The power measured at the output of an ideal band-pass filter fed by white noise is an 
estimate of the variance of the band-limited white noise (because the noise has zero mean 
value), and thus the power is a function of chi-square with n-l or 2BjT-l degrees of 
freedom. The filter power is assumed to be measured using an averaging time T. 

Therefore S(f c ) is given by 


2 

s<y = s<f c ) ^ <i7) 

The probability of S(f ) being within a certain interval of S(f ) can be calculated as fol- 

L c 

lows: 
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Let 


P[Si < S(f c ) < S 2 ] = l-a 


(18) 


where 


P[ a < y < b] 
1 - 0 ! 

S 1’ S 2 


probability of a < y < b 
desired confidence level 
to be determined 


From equations (17) and (18) 


2 

S 1 < s ( f c ) — < S 2 
n-1 


= 1-a 


Rn-l)S 

^) 


2 (n-l)S 2 

- x n-l - \ 


= 1 -a 


s(f c ) 


(19) 


Let 





1 - 0 ! 


( 20 ) 


where is the point of the i distribution that has a/2 area to the left and X 2 is 
the point that has a;/2 area to the right (X^ and X 2 can be found in tables of the chi- 
square distribution). Then, from equations (19) and (20) 


S l = 


S 2 " 


X l S ( f c> 

n-1 

X 2 S ( f C ) 

n-1 


( 21 ) 


If we now define the normalized power spectral density estimate as S(f c )/S(f c ), the fol 
lowing equation results: 


P 


n-1 “ S(f c ) “n-1 


1 - 0 ! 


( 22 ) 
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This equation, together with a table of the chi-square probability distribution, can be used 
to determine the normalized power spectral density at the tabulated confidence levels. 

As an example consider a = 20 percent 

B f = 10 Hz 

T = 1. 55 sec 

Then, the confidence level 1-a is 80 percent, and 

n-1 = 2B f T- 1 = 31-1 = 30 

2 

The values from v tables with a = 20 percent are 
(30) 

X 1 = 20. 60 
X 2 = 40. 26 


Thus 


P 


0. 686 < 


s(y 

s(f c ) 


< 1. 34 


= 80 percent 


For this numerical example, the probability of the measured spectral density lying in the 
interval [0. 686 S(f ), 1. 34 S(f )] is 80 percent. 

^ A 

It is unwieldy to use equation (22) to estimate the interval into which S(f c )/S(f c ) falls 
1-a percent of the time. A much simpler expression is obtained by using the normal ap- 
proximation to the chi-square distribution. For n-1 >30, the chi-square distribution is 
approximated very closely by the normal distribution with mean jli and standard devi- 
ation a: 


Let 


p = n-1 

a = y^(n- 1) 


J 


(23) 
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(24) 


s<g 

s(f c ) 


(n-l) 



where 

w is the normal random variable, with mean / a and standard deviation a. 

One can find in normal probability tables the values of 77 that satisfy 

P[- 7 y < z < 77 ] = 1 -a (25) 

where z is the standardized normal random variable (ju = 0 and a - 1). Now, 

w - a 
!± = z 

a 


and equation (25) becomes 


pU < — ~ ^ < 77 = 1 -a 

L a 


Using equations (23) and (24) in equation (26) results in 

r AT s(fj 

!-W < < 1 + fJlj - 

f n-l S(f c ) - 

or 



= 1 -a 


n-l 


-S(f c ) - 


= 1 -a 


where the percent error e is 


= = V J— - 

fn-1 f 2B f T- 1 


(26) 


(27) 


(28) 
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Figure 3. - Percent error as a function of the degrees of freedom for various 
confidence levels using the normal approximation to chi-square distribution 
(eq. (27)). 


For B f T » 1 

€ 88 ~ ~~ — (29) 

Equation (27) is to be interpreted to mean that the normalized values of the estimated 
spectral density will lie in the interval [(1 - e), (1 + e)] approximately l-a percent of 
the time, with e as given by equation (28). 

Equation (27) is plotted in figure 3. In the region to the left of (2BjT- 1) = 30, the 
graph gives only an approximation to the true confidence interval and is meant to be used 
as a general guide. For (2BjT-l) = 30, the interval limits are accurate to 1. 1 percent for 
a = 0. 50 and 11. 8 percent for a = 0. 05. If the normal approximation to the chi-square 
distribution is used in the previous numerical example, the following values are obtained: 

77 = 1 . 28 


and 


P 


0. 670 < 


S < f c> 

s(f c ) 


< 1.33 


= 80 percent 
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which agrees with the limits obtained previously within about 2. 3 percent. 

In summary, equation (22) is the exact expression used to obtain a confidence interval 
on the estimated spectral density when a confidence level is given. For large values of 
B f T, equations (27) and (28) may be used. 

The reader may recall equation (2) and might wonder why the response time of the 
analyzing filter did not enter into the analysis for random signals. The answer lies in 
the requirement that a number of independent samples of the power level be taken in each 
frequency interval. If the samples are to be truly independent, the filter output must 
have had ample time to change from one level to another. Thus, taking 2B|T-1 independ- 
ent samples in each frequency interval guarantees that the filter response time has been 
exceeded many times. 

Now that sufficient background information has been introduced the design of spectrum 
analyzers will be discussed. 


DESIGN 

The theory of swept frequency spectrum analyzers will be applied first to the design 
for periodic signals. Then the theory will be applied to random signals, which are rep- 
resentative of flight vibration data. 


Design for Periodic Signals 

For periodic signals, there is no statistical uncertainty in the time function and equa- 
tions (1), (2), and (3) apply: 


B f < ■ 


7TB, 


T d >> - 

f c 


( 1 ) 


( 2 ) 


( 3 ) 
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Equation (1) is satisfied by choosing a narrow band-pass filter as the analyzing filter. 
Equation (2) is satisfied by requiring each Fourier component to remain in the filter 



long enough to allow the filter to come up to nearly full response. Let the sweep rate of 
the analyzer be such that in time NjT c the filter moves across a frequency span of width 
B f (Nj^ is a numerical constant). If Nj = 3, the filter will be at approximately 95 percent 
of its full value as the Fourier component leaves the filter. For a constant sweep rate, 



(30) 


where 

f^ signal frequency being analyzed 

dfi 

— t sweep rate 

dt 

sweep rate constant 

Sweeping at a rate slow enough to allow adequate filter response and yet having minimum 
sweep time requires that 


Thus, 


N l T c 



dt 


B, 


K, 


N l T c 



(31) 


(32) 


Solving equation (30) subject to L (0) = f-^, where f ^ is the lowest frequency to be ana- 
lyzed, yields 

f^t) = Kjt + f L (33) 

The total sweep time t may be found by using L(t g ) = fjj, where fjj is the highest fre- 
quency to be analyzed, in equation (33). Then solving for t gives 

s 


15 


( 34 ) 


S K 1 ,B f 2 

Equations (32), (33), and (34) characterize the constant sweep rate spectrum analyzer for 
periodic signals. 

The restrictions imposed by equation (3) have not yet been incorporated into the 
design equations. Assume that the filter center frequency f in a swept frequency ana- 
lyzer with modulation is greater than 10 kilohertz. Then, from equation (3), the detector 
averaging time must be much greater than 100 microseconds. This criterion is easily 
met and in practical cases imposes no further restrictions on the sweep rate as given by 
equation (32). 

The design of a constant sweep rate spectrum analyzer for periodic signals is 
straightforward. A filter bandwidth Bj is chosen to give what the user considers to be 
adequate resolution, with equation (1) used as a guide. The sweep rate is given by equa- 
tion (32), the frequency-time function to be implemented by a variable frequency oscil- 
lator is given by equation (33), and the total sweep time is given by equation (34). 

2 

Notice in equation (32) that the sweep rate is directly proportional to B f , while the 

2 * 

total sweep time in equation (34) is inversely proportional to B^ . As a numerical ex- 
ample consider the following: 


B f = 10 Hz 
f R = 2000 Hz 
f L = 100 Hz 


Choose 


Then, 


Nj = n 


Kj = 100 Hz/sec 
f.(t) = lOOt + 100 
t g = 19 sec 
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This example has good resolution but a fairly long sweep time. If Bj were expanded to 
20 hertz, the total sweep time would be only 4. 8 seconds. The resolution would not have 
been degraded significantly; however, the system still would not meet the desired speci- 
fications listed in the INTRODUCTION because the minimum resolvable bandwidth would 
have been increased to 20 percent of the bandwidth center frequency at f 

Note that this analysis only applies to periodic data where there is no statistical un- 
certainty in the filter output. Spectrum analyzers for random data are discussed in the 
next section. 


Design for Random Signals 

Basic s weep rate equation . - Consider analyzing a spectral peak of width B n with a 
filter of width B^. First, to limit the amplitude error to e with a confidence level 1-a, 
the filter must be entirely within the peak B r for a time T determined from 

e=v + l—^— ( 28 ) 

y2B f T- 1 

The use of this equation assumes that 2B^T-1 is sufficiently large that the normal ap- 
proximation to the chi-square distribution may be used. 

Secondly, it can be shown that continuously sweeping a filter through a spectrum 
yields approximately the same results as moving the filter through it in discrete steps 
(ref. 7). From these considerations the filter can only be swept B n - Bj in frequency 
in a time T and still remain entirely within the peak for the duration T. Thus, 


gi , Af i . B n - B f 
dt At T 


(35) 


This is the basic sweep rate equation. 

Optimi zation of the filter bandwidth . - For in-flight vibration data the mechanical Q 
was assumed to be independent of the resonant frequency. Thus, if the minimum resolv- 
able peak with center frequency f^ is defined as B n , then the following equation results 
from the universal resonance curve: 

B = - A_ (36) 

^max 
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Clearly, if a filter of bandwidth is to resolve a spectral peak of width B n , then Bj 
must be less than B n . Let 

B f = C B n (37) 

where C is a constant to be determined. For a fixed bandwidth filter 


B. 


' C ( B n): 


Cf T 


min Q 


(38) 


max 


since the filter must be able to resolve a spectral peak at the frequency f^. Consider 
the optimum value for Bj or C. This means given € and 77 find B^ such that dL/dt 
is a maximum; that is, solve for 


d 

dB f 



= 0 


and 



< 0 


Solving equation (28) for T and substituting in the basic sweep rate equation yield 


Now 


Thus, 


df. 
1 

dt 


(B n - B f )2B f e' 

277 2 + € 2 


_d_A\ _ (B n -2B f )2 £ 2 
dB f \dt/ 277 2 + e 2 


(39) 
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B 

B f = — or C = - 
1 2 2 


Since 


_diri\ 4c 2 

d Bf 2 W a , 2 ** 2 


the sweep rate is truly a maximum. Note that C is a constant independent of e and 77. 

This optimization of B^/B n ignores blurring error. ^ However, consideration of the 
blurring error shows the apparent increase in bandwidth to be comparable to the resolu- 
tion and the frequency accuracy as given in the list of desired specifications (p. 3). For 
Bj/B n = 1/2 (the optimum ratio shown above), Ratz (ref. 6) showed that if a fairly sharp 
filter is used, the apparent increase in bandwidth is about 10 percent. 


Application 

Linear sweep (analyzer I) . - Using the basic sweep equation 


*1 _ B n - B f 
dt T 


(35) 


the desired types of sweeps can be generated. Assume B n , Bj, and T fixed. Thus, 


B - - B 
n 2 n 

dt T 


2Q T 
^max 


or 


f. 

1 


f L t 


2Q T 
^max 


+ 


(40) 


■*The blurring error is defined by Ratz (ref. 6) as an apparent increase in the band- 
width of the spectral peaks, the increase being due to the nonzero bandwidth of the filter 
and the shape of the transfer function. 
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This equation is for a linear sweep spectrum analyzer, which for brevity is called ana- 
lyzer I in this report. Such a system has constant frequency resolution rather than 
constant-percent frequency resolution; thus, when analyzing vibration data, much time is 
wasted scanning the high-frequency end of the spectrum. 

Nonlinear sweeps (analyze rs II and m). - If B n = f./Q mav as previously stated and 
Bj remains fixed, then the sweep rate equation becomes 


f i 


*1 

dt 


Q 


- 


max 


or 


f i B f^max + " B f^max^ e 


TQ 


max 


(41) 


In this case the frequency is an exponential function of time. Thus the system is nor- 
mally called an exponential sweep spectrum analyzer; however, it is called analyzer II 
in this report. 

As stated previously, this system is optimum for a fixed bandwidth filter and 
= f./Q ; however, the filter bandwidth is optimum only at f. = f T , where 
Bj = fe n ) min /2. Thus at the high end of the spectrum time is wasted because the filter is 
too narrow. This leads to the following considerations. Assume the existence of a vari- 
able bandwidth filter (ref. 14). Then the optimum filter bandwidth is 


1 2 2Q 


max 


Again the sweep rate equation is 


dfi 


dt 



T 


However, B f and T are no longer constants; hence, B^ and T must be substituted 

into the sweep rate equation in terms of f . . Using equation (28) and B f = f . /2Q 

i ii max 

yields 
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2,2 

€ f i 


Q^ax ^ 2 + **> 


or 



, 2 , 

f L e t 


Qmax ^ 2 + 2e2 > 


( 42 ) 


This system will be subsequently designated analyzer in. Note that the sweep rate is 
2 

proportional to T . For a comparison of the above systems consider the following ex- 
ample with 

B f = 5 Hz 
f L = 100 Hz 
f R = 2000 Hz 
e = 10 percent 

1-a = 90 percent 


Q 


max 


10 


From the normal distribution curve and the preceding value for 1-a, rj = 1. 64 in equa- 
tion (28); then substituting e = 10 percent into equation (28) yields T - 54. 2 seconds. 
Solving equations (40), (41), and (42) for the sweep time in all three systems yields the 
following equations: 


Analyzer I: 


l s * 


< E H- £ L> T < f H- f L)(^ 2+e2 ) 


B, 


2B 2 e 2 


(43) 
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Analyzer II: 


Analyzer III: 


t =Q 


max 


T In 


a H 


- B,Q 

f^max 

- B,Q 

f^max 


Qmax < 2 1 2 * f H - B lQmax 

2e^Bj f L ~ B f^max 

«. _ 2 < f H - f L)<ax (V + € 2 ) 


( 44 ) 


(45) 


Substituting the appropriate values of f H , f L , T, B f , Q max > V> 311(1 € into equations 
(43), (44), and (45) yields 


t g = 20 600 sec for analyzer I 

t g = 1980 sec for analyzer n 

t = 1030 sec for analyzer in 

s 

The nonlinear sweep systems have much shorter sweep times than the linear system; 
however, none are feasible for real-time analysis of in-flight vibration data with the ac- 
curacy and sweep time requirements given in the INTRODUCTION. Clearly, if t g is to 
be on the order of seconds, then accuracy, resolution, and confidence must be sacrificed. 
The following example was chosen to show to what extent the requirements must be de- 
graded to obtain sweep times on the order of seconds. Assume 

7 j = 0. 84 (i. e. ) 1-a = 60 percent 

f L = 100 Hz 

f H = 2000 Hz 


e = 40 percent 


22 



For analyzer I 


t = 91. 6 sec 


For analyzer II 


t = 8. 9 sec 
s 


For analyzer in 


t = 4. 6 sec 
s 

In the first numerical example, the value for the degrees of freedom is 540; so there 
is no doubt about the normal distribution being a good approximation to the chi-square dis- 
tribution. In the second example there are 9 degrees of freedom. The normal approxi- 
mation gives 


P 


0. 60 < 


S« c ) 
S(f c ) 


< 1.40 


0. 60 


The chi-square distribution gives 


0. 598 < 


S(f c > 
S(f c ) 


< 1.36 


0.60 


In view of the fact that e and a have been chosen to be 40 percent, the agreement be- 
tween the normal approximation and the chi-square distribution is more than adequate. 


BANDWIDTH REQUIRED FOR TRANSMISSION OF SPECTRUMS 
Single Spectrum Bandwidth Requirement 

Transmitting the spectrum of a signal rather than transmitting the time function will 
lead to a considerable saving in bandwidth. Consider an airborne spectrum analyzer with 
a resistor- capacitor (RC) averager on the power detector. Then, according to the dis- 
cussion following equation (13), RC = T/2. But according to equation (28) 
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( 46 ) 


,22 

T = 2t 7 + e 

2B f e 2 

For good response, the transmission channel should have a time constant shorter than the 
time constant of the data entering the channel; that is, 


data time constant 
channel time constant 


(47) 


The time constant of the channel response is 1/irB from equation (2), where B„ is the 

U C 

transmitting channel bandwidth. Substituting these results into equation (47) yields 


RC 

1/.B C 


= r 


or 

4rB f e 2 

B_ = i (48) 

L 9 9 

7t(2t 7 + e ) 

This expression is valid for all three types of sweep functions. For analyzers I and 
II, Bj is constant and equal to f T /2Q y , while for analyzer in the maximum value of 
B f is used in equation (48); that is, B f max = f jj/2Q max - As a result, analyzer HI re- 
quires substantially more bandwidth than analyzers I and II. The reason for the different 
bandwidth requirements is best understood by considering the power detector averaging 
time. In analyzers I and II, the averaging time of the power detector is constant, while 
in analyzer m the averaging time is inversely proportional to frequency. Thus, at the 
higher frequencies the power detector in analyzer HE responds much faster than the de- 
tector in the other systems. Providing the bandwidth to transmit the more rapid changes 
in power level from analyzer HI accounts for the increased bandwidth. The channel band- 
width B c is calculated for all three types of analyzers using the values of Bp e, and rj 
of the previous sweep time examples. The values for B c are given in table I for the 
following conditions: 
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TABLE I. - NUMERICAL EXAMPLES 


Analyzer 

Total sweep 
time, 

‘s’ 

sec 

Transmitting 

channel 

bandwidth, 

B c> 

Hz 

Telemetry 
bandwidth 
per channel, 

b m 

, 

m 

Hz /channel 

Telemetry 
bandwidth 
per channel 
for PCM system, 

M , 

■ \ m / PCM 
Hz/channel 



Example 1 


I 

20 600 

0.037 

0.082 

0.22 

n 

1 980 

.037 

.082 

.22 

m 

1 030 

.74 

. 164 

4.4 



Example 2 


i 

91.6 

i— i 

9. 1 

25 

ii 

8. 9 

4. 1 

9. 1 

25 

m 

4.6 

82 

182 

500 


Example 1 Example 2 


Lowest frequency analyzed, f ^ Hz 100 100 

Highest frequency analyzed, fjj, Hz 2000 2000 

Percent error, e 10 40 

Confidence level, 1-a, percent 90 60 

Maximum mechanical Q, Q 10 5 

Standardized random variable, rj 1.645 0.84 

Ratio of data time constant to channel time constant, r w ir 


Example 1 in the table is impractical because of the long sweep times. In example 2, B c 
for all the analyzers shows that spectrum transmission results in a substantial bandwidth 
saving over time function transmission; that is, 1900 hertz of data has been reduced to 
less than 100 hertz. 


Multispectrum Bandwidth Requirements 

Consider the transmission of the spectrums from many channels of wide-band data. 
Frequency multiplexing or time division multiplexing might be used to combine the spec- 
trums. Since the transmitting channel bandwidths, as calculated in the numerical ex- 
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amples, are so small, frequency stacking of the spectrum channels would be inefficient. 
More bandwidth would be used for channel separation than for information, and narrow, 
complex filters would be required. Time division multiplexing is more easily imple- 
mented. 

Consider the following sampling method. Assume each spectrum analyzer output to 
be band limited to frequencies less than B . Then from sampling theory, (ref. 12), the 
theoretical minimum sampling rate is 2B . Thus, to multiplex m spectrum channels, 

*” o c 

2mB samples per second are required. The amplitudes of the 2mB samples per sec- 
c c 

ond for m spectrum channels can be transmitted with an accuracy of about 3 percent; 

O Fj 

that is, e” * by using a telemetry channel whose bandwidth Bj^ is calculated from 

3. 5 (channel time constant) = ^ = — - — (49) 


or 


7mB, 


B M = 


(50) 


The quantity B^/m, the telemetry bandwidth per spectrum channel, may be used as a 
figure of merit for comparisons to other data transmission systems: 



m 



28 r B f € 
n 2 2rj 2 + e 2 


(51) 


Values of B M /m were calculated for the two examples, and the results are given in 
table I. These results show that even after multiplexing, a large bandwidth saving still 
remains. 


2 

The use of the theoretical minimum sampling rate based on the frequency limit B 
is sufficient in a practical system because of the frequency-limiting effects of the RC 
averager used on the power detector of the spectrum analyzer. The cut-off frequency of 
the RC averager is B c /jt. Thus, the frequency components in the analog signal repre- 
senting the spectrum are attenuated at frequencies above B c /tt. The lower frequency 
components of the signal, which comprise the major part of the signal, are sampled many 
times per cycle, which in a practical system yields good reproduction of the sampled 
spectrum. 
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If each sample were to be converted into a five- bit digital word, which would give an 
accuracy of about 3 percent, and an extra bit were added for synchronization, the bit rate 
for coding m spectrum channels would be 12mB c bits per second. It is a well-known 
rule in binary pulse transmission that through a low-pass channel that transmits up to a 
frequency of f d hertz, one can send 2f d pulses per second. Thus, for this binary- 
coded sampling scheme (pulse code modulation) the multiplexed bandwidth per spectrum 

channel [B 1 . /r /m ) is given by 
\ M /PCM 


*“] 


= 6B, 


m 


PCM 


24r B f e 2 
jt(2t 7 2 + e 2 ) 


(52) 


Values of (B.»/m) are also contained in table I. The bandwidth equations (48), (51), 
\ M /PCM 

and (52) take into account only the bandwidth required for the transmission of the spec- 
trums. Any low frequencies not included in the spectrum analysis must be transmitted 
separately. 

On the basis of these bandwidth calculations, analyzer II appears to be better suited 
for airborne use. In the two examples, analyzer n has a large bandwidth advantage 
(a factor of 20) and has a total sweep time less than twice that of analyzer m. Analyzer I 
has limited use because of its extremely long sweep times. 


CONCLUSIONS 

The optimum frequency- time functions were developed for three types of swept fre- 
quency spectrum analyzers. The functions were optimized by tailoring the design to the 
characteristics of the expected signals. The analysis for the linear and the exponential 
sweep analyzers (analyzers I and II, respectively) can be applied to systems that have 
already been developed. The frequency- time equation for the second-order sweep- rate 
spectrum analyzer (analyzer III) can be used as the basic design equation for the develop- 
ment of an improved spectrum analyzer having a short sweep time. 

The performance of in-flight spectrum analyzers on random data was calculated to be 
far short of the performance specifications desired. In order to obtain sweep times of the 
order of seconds, decreased resolution and/or reduced accuracy are required. 

The exponential sweep spectrum analyzer was judged to be the best spectrum ana- 
lyzer for in-flight use. Compared with the analyzer with the second-order sweep rate, 
the exponential sweep analyzer has only moderately longer sweep times (about two to one), 
requires less bandwidth (a ratio of 20 to 1), and is less complex. Both of the spectrum 
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analyzers with time-dependent sweep rates have sweep times about 10 times shorter than 
the constant sweep rate analyzer. 

Finally, the bandwidth per channel was calculated for the transmission of spectrums. 
Spectrum analysis showed substantial bandwidth savings over the transmission of wide- 
band data (at least 75 to 1). Many channels of spectrums can be transmitted over one 
wide-band channel of an IRIG FM/FM system. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 20, 1966, 

125-24-03-03-22. 



APPENDIX - SYMBOLS 


B transmitting channel bandwidth for one channel of continuous spectrum infor- 

c 

mation 



effective filter bandwidth (bandwidth of the ideal rectangular filter that passes 
the same power as the real filter when both are excited by white noise) 

telemetry bandwidth for m channels of sampled spectrums 
multiplexed bandwidth per spectrum channel for PCM system 


B n minimum resolvable spectrum peak with center frequency fj 

C ratio of B-. to B„ •„ 

I n, min 

E[y] statistical expected value of y 

2 

E[y ] statistical second moment of y 

fa’ fb lower and upper frequency limits on effective bandwidth of analyzing filter 

f analyzing filter center frequency 

fjj highest frequency to be analyzed 

f^ signal frequency being analyzed 

df 

— - sweep rate 

dt 

f ^ lowest frequency to be analyzed 

f mixing frequency 

s 

sweep rate constant 


m number of channels 

N number of samples taken to determine sample variance 

Nj numerical constant that allows adequate filter response for periodic signals 

n number of independent samples obtainable from band-limited white noise 
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P[ a < y < b] probability of a<y<b 

P ab average power in x(t) contained between frequencies f and f b 

P ab average power out of filter of bandwidth Bj, spanning f a to f b 

p[x(t)] probability distribution of x(t) 

Q mechanical Q 

R(t) autocorrelation function 

RC time constant of resistor-capacitor averaging circuit 

r ratio of data time constant to channel time constant 

S(f) power spectral density 

S(f ) exact value of power spectral density at f 

L v 

S(f c ) estimated value of power spectral density at f 

S ab (f) estimated value of average power spectral density between frequencies f & 

and f b 

S 2 lower and upper limits on confidence interval 

s^ sample variance 

T limit on integrations with respect to time; also used as averaging time for 

power detector 

T c time constant of filter response to sine wave 

T d detector averaging time 


T signal period for periodic signals 

ir 

t time 

t total sweep time 

Var x(t) variance of x(t) 

w normal random variable with mean p and standard deviation cr 

Xj, X 2 specific values of chi-square random variable; the values are those which 

give the desired confidence level 

x(t) signal (random or periodic) 
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Xj(t) filter output 

z standardized normal random variable with zero mean and a standard deviation of 

one 

a fraction of measured values of the power spectral density that will lie outside the 

confidence interval 

e percent error 

rj standardized random variable that gives desired confidence level 

li mean value 

2 

ct true variance 

t lead time in autocorrelation computation 

2 

chi-square random variable with N-l degrees of freedom 

Subscripts: 
max maximum 

min minimum 
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